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Abstract 

We give a hierarchy of many-parameter families of maps of the interval [0, 1] with an 
invariant measure and using the measure, we calculate Kolmogorov-Sinai entropy of 
these maps analytically. In contrary to the usual one-dimensional maps these maps 
do not possess period doubling or period-n-tupling cascade bifurcation to chaos, but 
they have single fixed point attractor at certain region of parameters space, where 
they bifurcate directly to chaos without having period-n-tupling scenario exactly at 
certain values of the parameters. 

1 Introduction 

In the past twenty years dynamical systems, particularly one dimensional iterative maps 
have attracted much attention and have become an important area of research activity. 
One of the landmarks in it was introduction of the concept of Sinai-Ruelle-Bowen (SRB) 
measure or natural invariant measure [Q, ||. This is, roughly speaking, a measure that 
is supported on an attractor and also describes statistics of long time behavior of the 
orbits for almost every initial condition in corresponding basin of attractor. This measure 
can be obtained by computing the fixed points of the so called Frobenius-Perron (FP) 
operator which can be viewed as a differential-integral operator. The exact determination 
of invariant measure of dynamical systems is rather a nontrivial task taking into account 
that the invariant measure of few dynamical systems such as one-parameter family of one- 
dimensional piecewise linear maps 0, || including Baker and tent maps or unimodal 
maps such as logistic map for certain values of its parameter can be derived analytically. 



Copyright © 2002 by M A Jafarizadeh and S Behnia 



Hierarchy of Chaotic Maps with an Invariant Measure 



27 



In most cases only numerical algorithms, as an example Ulam's method || Q, H, are used 
for computation of fixed points of FP-operator. 

Here in this paper we give a new hierarchy of many-parameter families of maps of the 
interval [0, 1] with an invariant measure. These maps are defined as ratios of polynomials 
where we have derived analytically their invariant measure for arbitrary values of the 
parameters. Using this measure, we have calculated analytically Kolomogorov-Sinai (KS) 
entropy of these maps. It is shown that they possess very peculiar property, that is, 
contrary to the usual maps, these do not possess period doubling or period-n-tupling 
cascade bifurcation to chaos, but instead they have single fixed point attr actor at certain 
region of parameters space where they bifurcate directly to chaos without having period- 
n-tupling scenario. 

The paper is organized as follows: In Section 2 we introduce the hierarchy of many- 
parameter families of chaotic maps. In Section 3 we show that the invariant measure is 
actually the eigenstate of the FP-operator corresponding to largest eigenvalue 1. Then 
in section IV using this measure we calculate KS-entropy of these maps for an arbitrary 
values of the parameters. Paper ends with a brief conclusion. 



2 Hierarchy of chaotic maps with an invariant measure 
and their compositions 

Let us first consider the one-parameter families of chaotic maps of the interval [0, 1] defined 
as the ratios of polynomials of degree N: 



$ N (x,a) 



a 



2 



( a 2 + {a 2 _ i)(_i)*r 2 Fi {-N, N, | , x) 
a 2 {T N {^)) 2 



(2.1) 



" l + (a?-l)(T N {y/xf) 
where N is an integer greater than 1. Also 

2 Fi (-N,N,~,x) = (-l) 7V cos(2iVarccos v / ^) = (-l) N T 2N (Vx) 

are hypergeometric polynomials of degree N and Tpj(x) are Chebyshev polynomials of 
type I || respectively. Obviously these maps unit interval [0,1] into itself. $n(x,o) is 
(N — l)-nodal map, that is it has (N — 1) critical points in unit interval [0, 1], since its 
derivative is proportional to derivative of hypergeometric polynomial 2-^1 (—N,N,^,x^ 
which is itself a hypergeometric polynomial of degree (N — 1), hence it has (N — 1) real 
roots in unit interval [0,1]. Defining Shwarzian derivative |i~0[| S^N(x,a) as: 



$%(x,a) 3 f^(x,a)\ 2 _ {<S>%(x,a)\ lf<f>%(x, a y 2 



& N (x,a) 2 \$' N (x,a)J \& N (x,a)J 2 \$' N (x,a 

with a prime denoting differentiation with respect to variable x, one can show that (see 
Appendix A): 

S ($jv(x, a)) = S Lf x (-N, N, ijj \ < 0. 
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Therefore, the maps $%(x) have at most N + 1 attracting periodic orbits [[h]]. Using the 
above hierarchy of family of one-parameter chaotic maps, we can generate new hierarchy 
of families of many-parameter chaotic maps with an invariant measure simply from the 

1,2,... ,n 



composition of these maps. Hence considering the functions (%, k 

<N 1} N 2 ,...,N n 



we denote their composition by: $ at,'^ "'°/v ; „ ( x ) wmc h can be written in terms of them in 



the following form: 



Ni,N2,...,N n 



(x) = ($jvi o $iv 2 o . . . o $jv n (ar)) 
^iVj • • ($iV„(iC,Q;n),a(n-l)) • • • , a 2 ), «i) ■ 



(2.2) 



Since these maps consist of composition of (iV& — l)-nodals (A; = 1,2, ... ,n) maps with 
negative Shwarzian derivative, they are (N1N2 ■ ■ ■ N n — l)-nodals maps and their Shwarzian 
derivative is negative, too [10|. Therefore these maps have at most N1N2 • • • N n + 1 
attracting periodic orbits [10|. As it is shown below in this section, these maps have only a 
single period one stable fixed points. Denoting m-composition of these functions by $( m ), 
it is straightforward to show that the derivative of <3>( m ) at its possible m x n periodic 
points of an m-cycle can be defined as: x^ k+ i = ®N k (x^k,ak), = $N n (x n ,fi., "at), 

and xi 1 = ®N n (xm,n,a n ), V = 1, 2 . . . , m, k = 1, 2, . . . , n is 



d_ 

dx 



$( m ) 



nn 

A*=l fc=l 



Nk 



[4 + (1 



a k) x V,k) 



(2.3) 



since for x^fc £ [0, 1] we have: 



min Q t + ( 1 — auX 



k x ^,k 



mm 



therefore: 



mm 



g?x 



I I — mm 



fc=l 



(1, 



<">7 



Hence, the above expression is definitely greater than 1 for Yi w~ < EI a fc < E[ ^fc> 

fc=i fc fc=i k=i 
that is, these maps do not have any kind of m-cycle or periodic orbits in the region of 

n n n 

the parameters space defined by Yl w~ ^ H a k < HI actually they are chaotic 

k=l k k=l k=l 

in this region of the parameters space. From (2.3) it follows that 



dx 



at 



m x n 



periodic points of m-cycle belonging to interval [0,1], vary between J\ (^k^k) 

k=l 



and 



k=i v 



for n a k < ri ih. and between ri (f^ 



k=l 



k=l 



k=l 



and J] ( N ka-k) m for ]J a k > 

k=l k=l 



Y[ N^, respectively. 
fe=i 

Definitely from the definition of these maps, we see that x = 1 and x = (in special 
case of odd integer values of N — 1, N2, ■ ■ ■ , N n ) belong to one of m-cycles. 



For n «* < rnuife 

cases in which x = 1 (x = 0) belongs to one of m-cycles, we have |-^-$( m )| < 1, hence 



II a k > n Nk ) , the formula (2.3) implies that for those 

k=l k=l 
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the curve of <]?( m ) starts at x = 1 (x = 0) beneath the bisector and then crosses it at the 
previous (next) periodic point with slope greater than one, since the formula (2.3) implies 
that the slope of fixed points increases with the decreasing (increasing) of therefore 
at all periodic points of n-cycles except for x = 1 {x = 0) the slope is greater than one 
that is they are unstable, this is possible only if x = 1 (x = 0) is the only period one fixed 
point of these maps. 

Hence, all m-cycles except for possible period one fixed points x = 1 and x = are 
unstable. 

Actually, the fixed point x = is the stable fixed point of these maps in the regions 

n n 

of the parameters spaces defined by > 0, k = 1,2, ... ,n and Y[ a k < Yl W~ om y f° r 

k=l k=l k 

odd integer values of N±, N 2 , ■ ■ ■ , N n , however, if one of the integers N^, k = 1,2, ... ,n 
happens to be even, then the x = will not be a stable fixed point anymore. But the 
fixed point x = 1 is stable fixed point of these maps in the regions of the parameters 

n n 

spaces defined by n a fc > Tl ^k an d a k < °Qj k = 1,2, ... ,n for all integer values of 

k=l k=l 

N 1: N 2 ,...,N n . 

As an example we give below some of these maps: 



a 2 (2x - 1) 



4x(l - x) +q 2 (2x - l) 2 ' 
a 2 x(4x - 3) 2 
a 2 x(4x - 3) 2 + (1 - x)(4x - l) 2 ' 



a 2 (l-8x(l-x)) z 



a 2 (l - 8x(l - x)) 2 + 16x(l - x)(l - 2x) 2 ' 



a 2 x(l6x 2 -20x + 5)' 



a 2 x(16x 2 - 20x + 5) 2 + (1 - x)(16x 2 - (2x - 1)) ' 
ai 2 (4x{x - 1) + (2x - l) 2 a 2 2 ) 2 
ai 2 (ix(x - 1) + (2x - l) 2 a 2 2 ) 2 - 16xa^(2x - l) 2 (x - 1) 
^ a2 {x) = ai 2 ((x - 1) (4x - l) 2 + x (4x - 3) 2 a 2 2 ) 2 
x (ai 2 ((x - 1) (4x - l) 2 + x(4x - 3) 2 a 2 2 ) 2 
-4xai(x- l)(4x- l) 2 (4x-3) 2 ) \ 



x ) = a 2 z [{x - 1) (4x - l) z + x (4x - 3) 2 q x 2 
x (u 2 2 ({x - 1) (4x - l) 2 + x(4x - 3) 2 ai 2 )' 
+ 4xa 2 (x- l)(4x- l) 2 (4x-3) 2 ) \ 



(2.4) 
(2.5) 
(2.6) 
(2.7) 

(2.8) 



(2.9) 



(2.10) 
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<f>^ Q * ( x ) = ai 2 a 2 2 x (4 x - 3) 2 ^3 (x - 1) (4 x - l) 2 + x (4 x - 3) 2 a 2 2 
- (x - l) 3 (Ax - l) 6 + 3x (3ai 2 - 2) {Ax - 3) 2 



x (x- l) 2 (4x- l) 4 a 2 2 + /i) , (2.11) 

where 

h = 3x 2 aj (-3 + 2a 2 ) (x - l)(4x - l) 2 (4x - 3) 4 + a 2 c^x 3 (4x - 3) 6 . 

Below we also introduce their conjugate or isomorphic maps which will be very useful 
in derivation of their invariant measure and calculation of their KS-entropy in the next 
section. Conjugacy means that the invertible map h(x) = maps / = [0, 1] into [0, oo) 
and transforms maps &N k (x,ak) into <b]y k (x,ak) defined as: 

&N k (x, ak) = h o & Nk (x, afe) o /i^ 1 ) = tan 2 (N^ arctan a/x) , 

a k 

Hence, this transforms the maps '"'°j^ n (x) ^jvViV2 "' £ JV„( :E ) defined as: 

;ai,a 2 ,...,a„ / s 
^N u N 2 ,...,N n V x ) 

= tan 2 ( Ni arctan ^/o—~ tan 2 ( N 2 arctan \/o ■ ■ ■ o — tan 2 (N n arctan ^/x) 
ai \ ai \ ai 



1 , 1 

—r tan iVi arctan 



07 



— ^ tan 2 ( A^2 arctan j / ■ ■ ■ — tan 2 (A n arctan y^) • • • ] I . (2.12) 
aS \ \ ai 'ii 



3 Invariant measure 

Dynamical systems, even apparently simple dynamical systems such as maps of an inter- 
val, can display a rich variety of different asymptotic behaviors. On measure theoretical 
level these types of behavior are described by SRB Jl], |llj or invariant measure describing 
statistically stationary states of the system. The probability measure /U on [0, 1] is called an 
SRB or invariant measure of the maps ^^^"'"n ( x ) §i ven m (2-2), if it is ^^'a^ "'"n ( x )~ 
invariant and absolutely continuous with respect to Lebesgue measure. For deterministic 
system such as these composed maps, the ^^'^''"'^ (x)-invariance means that its invari- 
ant measure u J) <«i,o2,...,a n (x) fulfills the following formal FP-integral equation: 

v N 1 ,N 2 ,...,N n 

^X".-.x (w) 1 6 (•" ^7;X.^S x)dx - 

This is equivalent to: 

^ Ql , Q2 ,..., Qn (y) = V] {j, ai , a2 ,..., an (x)^, (3.1) 

N 1 ,N 2 ,...,N n yy> 

defining the action of standard FP-operator for the map ^(x) over a function as: 

P s -i,« a -„ f(y) = J2 f( X )^ ( 3 - 2 ) 
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We see that, the invariant measure a d ,a 1 ,a 2 ,...,a n (x) is the eigenstate of the FP-operator 

^iVi ,JV,,...,JV.- 



P< 6 «i,a 2 >-.-,<*n corresponding to the largest eigenvalue 1. 

V N 1 ,N 2 ,...,N„ 



As we will prove below, the measure jx^i 



Afi ,N ,...,N n 



(x, (3) defined as: 



VP 



n y/x(l - x)(P + (1 - 0)x)' 



(3.3) 



is the invariant measure of the maps ^^'^''"''jv (x) provided that the parameter (3 is 
positive and fulfills the following relation: 



]Ja fc x 



k=i 



An - ^vZ^Zi 1 



A 



X • X 



V r >N2,N 3 ,...,N n {0 



2 V^.^S. -.^nUJ / 



where the polynomials A^ k {x) and B^ k {x) (k = 1,2, ... , n) are defined as: 



21 x ' 



1=0 



BN k {*)= E C 



2Z+1 ' 



«=0 



(3.4) 



(3.5) 
(3.6) 



where the symbol [ ] means greatest integer part. Also the functions 77^™ ( 4 J )^jv"_i'jVn ( ^ 1 j 



. . . and ^T/Vj'jvJ '"' jv„ ( /5 ) are defined in the following forms: 



C ( ^ 




a„_i,a„ 







a n -lA Nn _ 1 an 



a 2 ,a 3 ,...,a n j_ \ _ o 
r lN 2 ,N 3 ,...,N n \/3J~ P 



Oi2A.N 2 ( r)" 3 ' " 4 ' ■■■• an (±) 
1 \ r >N 3 ,N 4: ,...,Nn\t3) 



2 \ r lN 3 ,N 4 ,...,N n 



/ 



As we see the above measure is defined only for (3 > 0, hence from the relations (3.4), it 
follows that these maps are chaotic in the region of the parameters space which lead to 
positive solution of (3. Taking the limits of (3 — ► + and (3 — ► 00 in the relation (3.4) 

n n n 

respectively one can show that the chaotic regions is: ~f[ < Y\ a k < W ^k f° r odd 

k=l k k=l k=l 

integer values of N±, N2, ■ ■ ■ , N n and if one of the integers happens to become even, then 
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the chaotic region in the parameter space is defined by a k > 0, for k = 1,2, ... ,n and 

n n 

II a k < EI if one of the integers happens to become even, respectively. Out of these 

k=l k=l 

regions they have only period one stable fixed points. 

In order to prove that measure (3.3) satisfies equation (3.1), with (3 given by rela- 
tion (3.4), it is rather convenient to consider the conjugate map ^ < n^ < n^" ,C n„( x ) ^ w ^ n 
measure Uj«i,«2-,% denoted by us. related to the measure u a ,«i.«2- -.«n denoted by u$ 
through the following relation: 

Denoting ^ < N^ < N 2 " ,a N n ( x ) by V an< ^ inverting it, we get: 



2 / 1 / 9 ^l 71 " 

x kl = tan ( — arctan yj yaf + — 



' ' ! ''- air ' ( ~Jf 2 arctan \j x k x a\ + -J^j , 



Jt n 7T 



for kj = 1, . . . , A/j and j = 1, . . . , re. 

Then by taking the derivative of Xk 1>k2 ,...,k n with respect to y, we obtain: 



dx 



kl,kz,...,k n 



dy 



\k=i / 



Xk\,k2,--,kn 



y 



(1 + Sfcl,fc 2 ,-,fcJ( 1 + gfe 2 ,fc 3 ,-,fcn) •••(!+ gfc„-l,fcn)(l + SfcJ 

(1 + a&x fca>fcS) ... )fcf ,)(l + Q^-iXfcs,*!*,...,*:™) •••(! + ai^fe„_ 1 ,fe„)(l + a|x fcn )(l + afy) 



(3.7) 



to be substituted in equation (3.1). In derivation of above formula we have used chain 
rule property of the derivative of composite functions. 
Substituting the above result in equation (3.1), we have: 



fi*(v)Vv (1 + <*lv) = ( II tv| ) ' ' ' S y/ x kiM,..,kn 

\fc=l / fel fc 2 fc n 

(1 + SCfci,fc2,-,fcn)( 1 + X fc2,fc3,...,fcJ ■••(!+ ^fc„-l,fc„)(l + X k 



l + a 2 ) £*0«* 1) k 2) ... ) kj 



' (1 + alx k2t k 3 ,...,k n ){l + •••(! + ai^fe n _i,fc„)(l + a^fcj 

Now considering the following equation for the invariant measure /i$(y): 
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then the FP-equation reduced to: 



i,2 / n 

1 + a \V I TT a k 



n 



y, yv / (1 + Sfc 1 ,fc 2 ,...,fcJ(l + Sfc 2 ,fc 3 ,...,fcJ •••(!+ ^^fcJU + SfcJ 

J(l + ■(l + a§ a; fc n -i,0( 1 + Q! 2 a; fcn) 



Now using the follow identity: 



A? 



fc=0 f \0, 



for a one-parameter chaotic map y = &%{x) (Its proof is given in Appendix B), we obtain: 



JV« ^ 1 



k>2 >• • .fkn— 1 



C^i 1 + ^.te.--*- " ^ n (i) 1 + ^ U) x fel , fc2! ..., fcn _ 1 ' 



j9; i_r '/7VnV/3 

n "~ J ' n '' ^2 ^2 J ' k> ' A " A '" h ' r/ ' "" / ''" 1 



N n -iN n ^ f± (1 + a„-ia;fc 1 ,fc 2) ... ) fc f ,_ 1 )(l + 0x klM: ... u ) 
anan-i^„ (£) ^ (^Tj) i + a l_ lXkiM _ k _ 2 



(n \ N\ n2 

fc=i V fe 1= ife 2 =i 

/ (1 + Sfcl,fc 2 ,...,fcJ(l + gfc 2 ,fc 3 ,-,fcn) •••(!+ gfc„-l,fe w )(l + SfcJ 

a n-l X k3,k4,...,k„ 

) • • • (1 + a§x fcn _ life J(l + a^x fcn ) 



X 



k=i B Nn (±\ B N , 



n-l 



X • X 
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Now by inserting the right side of last relation in (3.5), we get: 

^ Ur.iXU)) B * («:k(ir)) 1 + (?) » 

We see that the above relation holds true provided that the parameter (3 fulfills the rela- 
tion (3.4). 



4 Kolmogrov— Sinai entropy 

KS-entropy or metric entropy [T1J measures how chaotic a dynamical system is and it 
is proportional to the rate at which information about the state of dynamical system is 
lost in the course of time or iteration. Therefore, it can also be defined as the average 
rate of information loss for a discrete measurable dynamical system (^^'^'"'^(x), fj). 
By introducing a partition a = A c (n\, . . . , n 7 ) of the interval [0, 1] into individual laps Ai, 
one can define the usual entropy associated with the partition by: 

71(7) 

Hfa-y) =-^m(i c )lnm(i c ), 

i=i 

where m(A c ) = j neA . [i(x)dx is the invariant measure of A{. By defining n-th refining 
7(71) of 7 as: 

n— 1 _fn 
fe=0 

then entropy per unit step of refining is defined by: 



1 



Now, if the size of individual laps of j(N) tends to zero as n increases, then the above 
entropy is known as KS-entropy, that is: 

KS-entropy, which is a quantitative measure of the rate of information loss with the refin- 
ing, may also be written as: 



dx ®N u N 2 ,...,N n W 



(4.1) 



which is also a statistical mechanical expression for the Lyapunov characteristic exponent, 
that is the mean divergence rate of two nearby orbits. The measurable dynamical system 
(^JVi'/Vb '"'"n ( a; )'A i ) * s chaotic for h > and predictive for h = 0. 



Hierarchy of Chaotic Maps with an Invariant Measure 



35 



In order to calculate the KS-entropy of the maps $^'#2 "' < tv„( x )' ^ ^ s rather convenient 
to consider their conjugate maps given by (2.12), since it can be shown that KS-entropy 
is a kind of topological invariant, that is, it is preserved under conjugacy map. Hence, we 
have: 



Using the integral (4.1), the KS-entropy of $ Ar 1 i '^'"' ^ n (x) can be written as: 



X 



o v/i(l + fix] 



In 



d 



tan 2 [N\ arctan JyN 2 ,N 3 ,-,N n ) 



( \ tan 2 (N 2 arctan y/yN 3 ,N 4 ,-,N n ) ■ ■ ■ -^-\tan 2 (N n arctan y/x) 

dyN :i ,N 4 ,;N n V«2 ' dxai 



or 



7T Jo 



+ 



o v^(l + ^) 

ffidx 



i r y^ 

7o V^C 1 - 



d 



— tan 2 (iVi arctan y/yN 2 ,N 3 ,;N n ) 
dyN 2 ,N 3 ,-,N n \a{ 



d 



+ ••• + 

+ 



+ /3x 
00 y^dx 



dyN 3 ,N4„-,N n \a> 



2 tan 2 (AT 2 arctan v / yjv 3 ,JV4,-,JV„) 



TT J s/x(l+fix) 
yfj3dx 



i r c 



v^(l + /3x) 



In 
In 



d 



1 



dVN n V«n-1 

d / 1 



tan 2 (JV n _i arctan y 7 ?/^) 



tan 2 (iV n arctan (\/x ) ) 



(4.2) 



where 



tan 2 (7V n arctan (a/x)) , 



VN„- U N n = tan 2 (iV n _i arctan( y / ^))) , 



a n -i 



(4.3) 
(4.4) 



yN 2 ,N 3 ,-,N n = — tan 2 (iVi arctan( v / yjv 3i jv4,-,JV n )) 



a 



(4.5) 



Now, we calculate the integrals appearing above in the expression for the entropy sepa- 
rately. The last integral is calculated in appendix C which reads: 



1 f c 
n Jo 



s/fid. 



X 



= ln 



In 



d ( 1 



dx \ a 2 , 



tan 2 [N n arctan (y^c )) 



An„ (j) £jv„ 



(4.6) 



36 



M A Jafarizadeh and S Behnia 



In order to calculate the integral before the last one in (4.2), that is the following integral, 



If 

71" JO 



\ff3dx 



y/x{l + fix 



In 



d 



tan 2 (N n -i arctan y/yN n ) 



(4.7) 



first we make the following change of variable by inverting the relation (4.3) 



f ' ,m ~ ( iJ— arctan (\ 7Zv~o 2 , ) + k " 



iVn-1 



JVn-1 



fc n _i = 1, . . . , iV, 



71- 1- 



Then the integral (4.7) is reduced to: 



J- 1 ^ «*- I. 



™n — 1 — 1 ^n— 1 



In 



1 



<&JV„ V"n-1 



tan 2 (7V n _i arctan y^y^) 



where x kn 1 and x£ (^n-i = 1,2, . . . , iV„_i) denote the initial and end points of fc-th 

branch of the inversion of function = tan 2 (JV n arctan y/x)^ respectively. Now, by 

inserting the derivative of Xk n _ 1 with respect to yjy„ in the above relation and changing 
the order of sum and integration, we get: 



^ -/o k ^ =l N n _ ly /y^{l + a z n _ iyNn ) 



^/x^I7(l + /3x fcn _ 1 ) 



x In 



d VN n \«n-l 



tan 2 (iV„_i arctan 



By using the formula (B.6) of Appendix B, we get: 



i r 

n Jo 



Vf3dy r , 



B N n ( \ 



a n A- 1 



g I U n S±N n I g 



x In 



d 



d VN n V«n-1 



tan 2 (iV n _i arctan ^/yv n ) 



Finally, through calculating the above integral with the prescription of Appendix C, we 
obtain: 

/ -v (l _j ( i + +2j//;r : " ! 

In 



te) ^ te) 



Similarly, we can calculate the other integrals appearing in the expression for the entropy 
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of the composed maps given in (4.2): 



7T J y/x{l + (3x) 



111 



h 
d 



N u N 2 ,...,N n 



1 



dVN k ,N k+1 ,;N n 



a 



In 



2 — tan 2 (iVfc_i arctan y/yN h ,N k+1 ,;N„) 

(JV fc _!)-l\ 



fc-1 



•,0!n / 1 

..,AT„ 1/3 



V 



^ (J)) ^ (* 



for k = 1,2, ... ,n. 

Finally, summing the above integral, we get the following expression for the entropy of 
these maps: 



In 



x + \ hZ 



2(JV„_i-l) 



{N l N 2 ---N n )[l + y//3 

_ 2(ATi- l ! 



2(JV„-1) 



(4.8) 



Using the formulas (4.8), one can show that KS-entropy of one parameter families has the 
following asymptotic behavior: 



h ((J,,<f>N {x,a = N + 0")) ~ (N- a)i 



h [ [J,, $ N [ x, a = — + (T 



~ a 



1 \ ' 

N 



near the bifurcation points, that is (3 — ► as (3 — ► oo. The above asymptotic behaviors 
indicate that one-parameter maps $%(x) belong to the same universality class which are 
different from the universality class of pitch fork bifurcating maps, but their asymptotic 
behavior is similar to class of intermittent maps [Oj. Even though, intermittency can 
not occur in these maps for any values of parameter a, the maps <&^(x) and their n- 
composition <l?( n ) do not have minimum values other than zero and maximum values other 
than 1 in the interval [0, 1]. 

Also by imposing the relations between the parameters a^, k = 1,2, ... ,n which are 
consistent with the relation (3.4), we can reduce these maps to other many-parameter 
families of maps with the number of the parameters less than n. Particularly by imposing 
enough relations, we can reduce them to one-parameter families of chaotic maps with an 
arbitrary asymptotic behavior as the parameter takes the limiting values. Hence we can 
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construct chaotic maps with arbitrary universality class. As an illustration, we consider 
the chaotic map ^^(x). Using the formula (4.8), we have: 

h (u *™(x)\ ~ ln (1 + ^ )2(2 ^ + " 2(1 + ^ 2 

with the following relation among the parameters a±, «2 and (3: 

ai (4/? + «i(l + /?) 2 ) =4q 2 /3(1 + /3) 

which is obtained from the relation (3.4). Now choosing f3 = a^, < v < 2, the above 
relation reduces to: 

4o^(l + a|) 

^ «2 (1+a , )2+4a , 

and entropy given by (4.8) reads: 

/ \ fl + af) 2 f2a| +Q 2 (l + a^) N2 

r' 2 ' 2V V (l + a^)(4a^ + a 2 (l + ^)2) 

which has the following asymptotic behavior near a 2 — ► and a 2 — ► oo: 

/i ^u, $2,2 ( x )) ~ a 2 as a 2 — ► 0, 
/i ^/i, <I>2 2 2 (2;)^ ~ (~~) aS 02 — * °°' 

The above asymptotic behaviours indicate that for an arbitrary value of < v < 2 
the maps $2 \{ x ) belong to the universality class which is different from the universality 
class of one-parameter chaotic maps of $n(x) (2.1) or the universality class of pitch fork 
bifurcating maps. 



5 Conclusion 

We have given hierarchy of exactly solvable many-parameter families of one-dimensional 
chaotic maps with an invariant measure, that is measurable dynamical system with an 
interesting property of being either chaotic or having stable fixed point, and they bifurcate 
from a stable single periodic state to chaotic one and vice- versa without having usual period 
doubling or period-n-tupling scenario. 

Again this interesting property is due to the existence of invariant measure for a region 
of the parameters space of these maps. Hence, to support this conjecture, it would be 
interesting to find the other measurable one-parameter maps, specially coupled or higher 
dimensional maps, which are under investigation. 



Hierarchy of Chaotic Maps with an Invariant Measure 



39 



Appendix 

A Shwartzian derivative 



The Shwarzian derivative S&n(x) [1C] is defined as: 



&»(x) 3 f$%(x) 



$%(x)\' 1 



with a prime denoting a single differential. One can show that: 



N-l 



S ($ N (x)) = S ^ 2 Fx \-N, N,-,x] ) < 0. 
since ^ (2-F1 [~N,N, \,%)) can be written as: 

i(^(-^,w,i,«))-Vn(«-^) 

with < xi < X2 < X3 < • • • < a?jv_i < 1, then we have: 

• 1 1 



5 I 2 Fi [ -N,N,l,x 



2'" J J 2 ^ {x- Xj f (x-xj) 



< 0. 



Also, one can show the shwartzian derivative of composition of function with negative 
shwartzian derivatives is negative too. 

B Derivation of the formula (3.9) 



In order to drive formula (3.9), we write the summation in its left side as: 

1 + c z 

N 



a ^ 1 + or Xk a ( (3 — 1 \ d 



(^)^HjrH) 



(B.l) 



To evaluate the second term in the right side of (B.l), we denote &%(x) by y therefore, 
the map (x,a) = -\ t&n 2 (N arctan -y/x) can be written as: 

= a 2 y cos 2 (iV arctan \fx) — sin 2 (iV arctan \fx). (B-2) 

Now, we can write the equation (B.l) in the following form: 

/ (m 



N 



1 + a z x k 



•l) i 



o 



[V] 



X 



fe=0 




N x k 



N 



Y[(x - X k ), 



fe=l 
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where Xk are the roots of equation (B.2) and they are given by 

(X k7r\ 
— arctan y ya 2 + — J , k = 1, . . . , N. 

Therefore, we have: 

d 



(B.3) 



80- 



-in (n^~ i+ ^)) 



d 



d/3- 



j In (1 - /T 1 )^ (a 2 y cos 2 (n arctan y 7 -/?" 1 ) - sin 2 (iV arctan y 7 -/?" 1 )) 



iV^_ , (1 + » 2 ^) -4 (fr) 



/3-1 



2 ' 



(B.4) 



with polynomials A(x) and -B(x) given in (3.5) and (3.6) where in derivation of above 
formulas we have used the following identities: 



cos {N arctan \/x) 



A(-x) 

N ' 

(1 +x)t 



sin (iV arctan \/x) 



B(-x) 
(1 +x)~ 



By inserting the results (B.4) in (B.l), we get: 



v 



a 1 + a 2 Xk 
l + px k 

fe=0 



1 + a 2 y 



B I 



«4 (A 



P 



aA 



B I 



(B.5) 



(B.6) 



C Derivation of the last integral of (4.2) 

Using the relations (3.5) and (3.6) the last integral of (4.2) can be written as: 

(l + x) JV " 1 J B(-x) 



* JO 



y/pdx , ( N 
m 



y/x(l + px) \a 2 



(A(-aO)S 



(C.l) 



We see that polynomials appearing in the numerator (denominator) of integrand appearing 
on the right side of equation (C.l) have ^ 2 ^ (^f) simple roots denoted by , k = 

1, . . . , [^f^] (a^fc, fe = 1, . . . , [y]) in the interval [0, oo). Hence, we can write the above 
formula in the following form: 



h(a ^Uir ^ In 



fJV-ll 



iV 



(i+x)^- 1 n k-^fc 1 



fe=i 



n N-^i 3 
fe=i 
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Now, making the following change of variable x = tan 2 ® , and taking into account that 
degree of numerators and denominator are equal for both even and odd values of N, we 
get: 



-| roc ) /AT 



+ (N - 1) In |/3 + 1 + (J3 - 1) cos 6\ 



+ J2 ln|l-xf/3+(l + xf/3)cos0|-3^1n|l-x^/?+(l + x^/?)cos0| \. 



k=l 



k=l 



Using the following integrals: 

i r 

— In \a + b cos 0| = < 
71" Jo 



In 



In 



a + V a 2 — b 2 



we get: 



h 



In 



|a| > \b\, 



N(l + [3 + 2^p) 



N-l 



\N-1] 



(e^I IE'*/ 

y y fc=o / \ fc=o 



(C.2) 
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